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Abstract 

Quantum dynamical semigroups are applied to the study of the time 
evolution of harmonic oscillators, both bosonic and fermionic. Explicit 
expressions for the density matrices describing the states of these sys- 
tems are derived using the holomorphic representation. Bosonic and 
fermionic degrees of freedom are then put together to form a supersym- 
metric oscillator; the conditions that assure supersymmetry invariance 
of the corresponding dynamical equations are explicitly derived. 
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1. INTRODUCTION 



The dynamics of a small system S in interaction with a large environment £ is in 
general very complex and can not be described in terms of evolution equations that are 
local in time. Possible initial correlations and the continuous exchange of energy as well 
as entropy between the S and £ produce phenomena of irreversibility and dissipation. 

Nevertheless, there are instances for which simple and mathematically precise descrip- 
tion of the subdynamics can actually be given. When the typical time scale in the evolution 
of the subsystem S is much larger than the characteristic time correlations in the envi- 
ronment, one expects (and actually proves) the disappearance of memory and non-linear 
phenomena, although quantum coherence is usually lost. [1-5] 

In such cases, the states of S, conveniently described by a density matrix p, are seen to 
evolve in time by means of a family of linear maps that obey very basic physical require- 
ments, like forward in time composition (semigroup property) and complete positivity. 
They form a so-called quantum dynamical semigroup. [1-3] 

This description of the time evolution of open quantum systems is actually very gen- 
eral; it is applicable to all physical situations for which the interaction between S and £ can 
be considered to be weak and for times for which non-linear disturbances due to possible 
initial correlations have disappeared. In particular, quantum dynamical semigroups have 
been used to model laser dynamics in quantum optics, [6-8] to study the evolution of vari- 
ous statistical systems, [1-3] to analyze the interaction of a microsystem with a macroscopic 
apparatus. [9-11] 

Recently, they have been used to describe effects leading to irreversibility and dissi- 
pation in elementary particle physics phenomena. Non-standard low energy effects accom- 
panied by loss of quantum coherence are in fact expected to appear as a consequence of 
gravitational quantum fluctuations at Planck's scale. [12] Detailed analysis of these effects 
have been performed in the system of neutral mesons, [13-16] in neutron interfere metry, [17] 
neutrino oscillations, [18] and in the propagation of polarized photons; [19, 20] the outcome 
of these investigations is that present and future elementary particle experiments will likely 
put stringent bounds on these non-standard dissipative phenomena. 

These studies, in particular those dealing with correlated neutral mesons, [14, 21] have 
also further clarified the importance of the condition of complete positivity in the de- 
scription of open quantum systems. In many investigations complete positivity is often 
replaced by the milder condition of simple positivity; this guarantees the positivity of the 
eigenvalues of the density matrix of the subsystem <S, but not that of a more general sys- 
tem obtained by trivially coupling S with an arbitrary n-level system. Lack of imposing 
this more stringent requirement could lead to unacceptable physical consequences, like the 
appearance of negative probabilities. [21] 

To further analyze the properties of the quantum dynamical semigroup description 
of open systems, we shall apply this general framework to the analysis of the evolution 
of one dimensional oscillators, both bosonic and fermionic (for earlier investigations on 
the bosonic case, see [22-25] and references therein). We shall adopt the holomorphic 
representation [26, 27, 23] since it allows an explicit description of the relevant density 
matrices in terms of complex and Grassmannian (anticommuting) variables; in the most 
simple situations, the general form of these density matrices turns out to be Gaussian. This 
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allows the explicit evaluation of the corresponding (von Neumann) entropy and analysis of 
its time evolution. Finally, we shall combine a bosonic and a fermionic degrees of freedom 
to form a supersymmetric oscillator. We shall then derive the conditions that guarantee 
the supersymmetric invariance of the dynamical equations and discuss how these affect the 
time evolution of the total density matrix. 



2. THE BOSONIC OSCILLATOR 

As explained in the introductory remarks, we shall study the dynamics of a single 
oscillator in interaction with a large environment. The states of the system will be repre- 
sented by a density matrix ps, i-e. by a positive hermitian operator, with constant trace, 
acting on the bosonic Hilbert space Hb- Our analysis is based on the assumption that its 
time evolution is given by a quantum dynamical semigroup; this is a completely positive, 
trace preserving, one parameter (=time) family of linear maps, acting on the set {p} of 
bosonic density matrices. These maps are generated by equations of the following general 
form: [1-3] 

dp(t) 



dt 



= jC[p(t)] = -i[H,p(t)}+L[p(t)} . (2.1) 



The first term in C is the standard quantum mechanical one, that contains the system 
hamiltonian H, driving the time evolution in absence of the environment. In the case 
of the bosonic oscillator, it can be taken to have the most general quadratic form in the 
bosonic creation and annihilation a operators: 

H B = - [u B (a 1 a + aa 1 ) + pa 2 + p* a) 2 ] , (2.2) 

where ujb > and p is a complex parameter (the star means complex conjugation). The 
second piece L[p] takes into account the interaction with the environment; it is a linear 
map, whose form is fully determined by the requirement of complete positivity and trace 
conservation: 

L[P] = -IJ2( L 1 L ^P + P L 1 L ^ +E L ^4- (2-3) 

k k 

The operators L k should be chosen such that the expression in (2.3) is well defined. In 
absence of the term L[p], pure states would be transformed into pure states. Instead, the 
additional piece (2.3) produce in general dissipation and loss of quantum coherence. 

The choice of the operators L k is largely arbitrary. However, since the hamiltonian 
Hb is quadratic in a* and a, one is led to assume the same property also for the additional 
term L[p\. This implies a linear expression for the operators L k : 

L k = r k a J r s k a) , (2.4) 

with r k and s k complex parameters; this requirement further guarantees the exact solvabil- 
ity of the equation in (2.1). Note that the operators (2.4) are not bounded; nevertheless, 
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by adapting the arguments presented in Ref.[22] to the present case, one can show that 
the exponential map generated by (2.1) is well defined. 

This description of the damped bosonic oscillator can be further simplified by means 
of a suitable canonical transformation. First, notice that not all values of the oscillator 
frequency u B and the complex parameter \x are physically allowed. Indeed, the spectrum 
of the hamiltonian in (2.2) is bounded from below only for: 

w|-|H 2 >0. (2.5) 

This is a consequence of the fact that H is an element of the Lie algebra su{l, 1), whose 
generators in the so-called metaplectic representation take the form: 

Ko = I( a taW) , (K+) t = K_ = y. (2.6) 

The condition (2.5) guarantees that H can be unitarily "rotated" to an element of the 
Cartan algebra with spectrum bounded from below. [23] 

In other terms, by means of a unitary canonical transformation, one can now pass to 
new operators: [26, 28] 

a = $a + ^a t , 
a* = ^* a + $* a 1 , 

with 



(2.7) 



such that the hamiltonian take the simplified form: 

H B = ^-{a\~a}. (2.9) 

The operators in (2.3) are still linear in the new variables a* and a, although with 
redefined coefficients. 

This discussion explicitly shows that, without loss of generality, one can set \i = in 
(2.2); a non vanishing \i can always be reinstated at the end by undoing the transformation 
(2.7). With this choice, the evolution equation (2.1) for the bosonic oscillator becomes: 

d -P^L = C B [p B {t)\ = -tu B [a^a lPB (t)}+L B [p B (t)} , (2.10) 
where, by inserting (2.4) into (2.3), one has: 



l b[p] =r] B [\ [ ap,a ] ] + [a^pa^ +<j B ([o t p,a] + [a f , pa]j 

- A B [a, [a, p\] - \ B [a f , [a 1 , p]] 



(2.11) 
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with 

VB = 2^l rfc ! 2 ' a B = -^\Sk\ 2 , ^B = ~J2 r k S k- (2-12) 

k k k 

Note that from these expressions one deduces that: 

VB>0, a B >0, \\b\ 2 <Vbvb, (2.13) 

the last relation being a consequence of the Schwartz inequality; let us remark that these 
are precisely the conditions that assure complete positivity of the time evolution generated 
by the operator Lb in (2.11). 

In order to study the solutions of the equation (2.10), we shall work in the holomorphic 
representation; [26-28, 23] it allows deriving explicit expressions for the density matrix ps{t) 
so that its behaviour in various regimes can be more easily discussed. In this formulation, 
the elements of the bosonic Hilbert space Hb are represented by holomorphic functions 
ip(z) of the complex variable z, with inner product:^ 

= J ^{ z )4>{z) e~~ zz dzdz . (2.14) 

To every operator O acting on Tis there correspond a kernel 0(z,z) of two independent 
complex variables z and z, such that for the state \<f>) = O \ ip) one finds the representation: 

cf)(z) = J 0{z,w)ip{w) e~™ w dwdw . (2.15) 

In particular, the creation and annihilation operators, when acting on a state |^), are 
realized by multiplication and differentiation by the variable z: 

aV>-zV(«), a|V>>^JU(*>, (2-16) 

while the identity operator is represented by e zz . 

Since the term Cb in (2.10) is at most quadratic in a} and a, the kernel pb{z, z;t) 
representing the solution of (2.10) can be taken to be of generic Gaussian form: 

p B (z,z;t) = _^ =e -^ m by^zz-x{t)z^- x {t)^]+zz (2 _ 17) 

Trace conservation for all times, 

Tr[p B (t)} = J p B (z,z;t)e- zz dzdz = \ , (2.18) 

readily implies: 

N{t)=y 2 {t)-\x{t)\\ (2.19) 



^ Here and in the following we use the conventions of Ref. [26] 



5 



while using (2.16) one finds that the unknown functions x(t), x(t) = [x(t)]* and y(t) have 
the following physical meaning: 



/q2 
— [p B (z,z;t)] e~~ zz dzdz = x(t) , 

(a t2 )(t) = Tr[a t2 p B (t)] = J z 2 p B (z, z;t) e~~ zz dz dz = x(t) , (2.20) 

/d 
— [z pb(z, z; £)] e~ zz dzdz = y(t) . 

For simplicity, in writing (2.17) we have assumed (a^)(t) = (a)(t) = for all times. As 
shown in the Appendix, this condition can be easily released starting with a more general 
Ansatz for pb(z, z;t); it will not be needed for the considerations that follow. 

Inserting (2.17) in the evolution equation (2.10), with the help of the relations (2.16) 
one finds that the unknown functions x(t) and y(t) satisfy the following linear equations: 

x(t) = -2 (t] B -<t b + iojb) x{t) - 2X B , 
y(t) = -2 ( VB - o B ) y(t) + 2t] B . 

General solutions can be easily obtained. For initial values xq = x(0), yo = y(0) and 
Vb 7^ cr_B, one finds: 

x(t)=E(t)e- 2iaJBt (x -x oo )+x oo , 
y(t) = E(t) (y - yoo) + y^ , 



(2.22) 



where 

E(t) _ e -2(vB-* B )t x _ *B(<rB-VB+iu> B ) _ VB 

{VB-o-Br + ^B Vb-cfb 

while in the particular case t\b = o~b- 

x(t) = e ~ 2luJBt (x - x c ) + x c , x c = i— , 



(2.23) 



ojb (2.24) 



y(t) = 2r] B t + y . 



The large time behaviour of these solutions depends on the relative magnitude of the 
two positive parameters r\B and as- Only when t]b > ob, the functions x(t) and y(t) have 
a well-defined limit. In this case, independently from the initial conditions, the density 
matrix p B (t) approaches for large t the equilibrium state , obtained substituting in 
(2.17) the asymptotic values Xqo and yoo for x(t) and y(t). Indeed, is clearly a fixed 
point of the evolution equation (2.10), for any value of t\b and as- 

Notice that does not correspond in general to a thermal equilibrium state; to 
obtain an asymptotic Gibbs distribution, one has to set Xb = and introduce the inverse 
temperature (3 via the condition y^ = [coth(f3uB /2) + l]/2 (compare with (2.27) and 
(2.28) below), or equivalently e^ B = 7] B /<Jb- 
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On the other hand, when t]b < o~b, the exponential term E(t) in (2.22) blows up for 
large times, while in the special case t/b = &b, y(t) grows linearly in time and x(t) has an 
oscillatory behaviour. In both cases, for generic initial conditions, the normalization factor 
N(t) in (2.19) grows unbounded, so that the functional pB(z,z;t) becomes vanishingly 
small, while retaining its normalization, Tr \pB(t)] = 

This peculiar behaviour can also be analyzed with the help of the Weyl operators: 



W\u] 



va-\-ua' 



(2.25) 



By studying the time evolution of W induced by (2.10) via the relation Tr[W(t) ps] = 
Tr[W ps(t)], one finds that when t]b > o~b all Weyl operators remain well-defined for all 
t, approaching the identity for large times; on the other hand, for t\b < o~b one discovers 
that all Weyl operators vanish in the large time limit, except the identity W[0], which is 
clearly a fixed point of the time evolution. 

As mentioned in the Introduction, the entropy of an open system usually varies with 
time, due to the interaction with the environment. In many physical instances, monotonic 
increase of the von Neumann entropy, 



%] = -Tr[p lnp] = -(lnp) , 



(2.26) 



is a desirable property. [1-3, 13-20] For the damped oscillator described by the evolution 
equation (2.10), this request can not be fulfilled in general, if one insists on the existence 
of a well-behaved large-time equilibrium limit. 

The explicit evaluation of Sb = S[pb] is simplified by noticing that to the represen- 
tation kernel (2.17) there corresponds the operator Ansatz: 



Pb = 



1/2 



[A(aa f +a f a)+Ba 2 +Ba 1 2 ] 



_coth 2 (fi/2) - 1 

where the parameters A and B are related to x and y of (2.17) through the relations: 



(2.27) 



x 



B , O 

-20 C ° th 2 



V= 2 



A , O 
C ° th 2 +1 



0= {A 2 -\B\ 2 ) 



2\V2 



Inserting these relations in the definition (2.26), one obtains: 



Sb 



-I In 
2 



.coth 2 (fi/2) - 1 



+ 2 



A (aa f + a f a) + B (a 2 ) + B (a t2 ) 



(2.28) 



(2.29) 



It is now convenient to introduce the following quantity: 



X 2 =-(aa^ + a^a) 2 - 



(a 2 ) (a 



t2\ 



< 2 /-l/2) 2 -|x| 2 = ^coth 2 ^ , 



(2.30) 



t For i]b < &b, the exception is given by = since, as noted before, Xoo and j/oo 

are fixed points of (2.21). Note that yoo blows up for t]b = &b, so that also becomes 
vanishingly small in this limit. As a consequence, the infinite temperature limit is singular. 



7 



such that % > 1/2; in terms of this variable, one can easily rewrite (2.29) as: 

S»=(x + i)ln(x+i)-(x-i)ln(x-i). (2.31) 

The entropy Sb always grows with x, starting at the minimum Sb = for x = 1/2 
and increasing as lnx for large x- Recalling the explicit time-dependence of x and y 
in (2.22) and (2.24), one realizes that in general for t\b ^ <?b the variable x does not 
monotonically grows with t, so that the condition Sb > can not be satisfied for all times. 
In particular, when t\b > &b the equilibrium state is reached in general at the expense 
of some negative entropy-exchange with the environment. 

The case t\b = &b is again special; in fact, the operators in (2.4) are now hermitian 
and therefore the condition 5^ > is guaranteed. Indeed, in this case x approaches infinity 
for large times, and therefore so does Sb- Alternatively, using (2.3) in the definition (2.26), 
one can directly show that: [29] 

Sb > ( J2 4] ) = 2(^ B - V b) ■ (2.32) 

k 

Note however that Lk = L\ is only a sufficient condition for entropy increase. Indeed, 
for t]b > ob take x$ = x^ and yo < y^; in this case x grows with time since y{t) is always 
positive, and therefore also Sb never decreases. 



3. THE FERMIONIC OSCILLATOR 

We shall now extend the analysis of the previous section to the case of a fermionic 
oscillator. The corresponding creation cr an( j annihilation a operators obey now the 
algebra: 

{a,a f } = 1 , a 2 = cJ 2 = . (3.1) 

As in the bosonic case, we shall assume the system in interaction with a large environment, 
and describe its time evolution by means of a quantum dynamical semigroup. 

The states of the system will be described by an appropriate density matrix pp, acting 
on the elements of the fermionic Hilbert space TCf- This operator obeys an evolution 
equation of the form (2.1), where now the hamiltonian can be taken to be: 

H F = ^- [a+, a] , u F > . (3.2) 

Since and a are now nilpotent, the additional piece L[p\ in (2.3) turns out to be at most 
quadratic in these variables, and the operators Lk assume the generic form 

L k = r' k a + 4 a ] . (3.3) 
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Inserting this in (2.3), one explicitly finds: 



L F [p] = rjp (2 apex* — a^ap — pa^a) + a F (2 pa — aa^ p — paa^) ^ 

+ 2(X F apa + X F a} pa^ , 

where the parameters rj F , &f and X F are as in (2.12) with the coefficients and Sk replaced 
by the primed ones. Then, the complete evolution equation for the density matrix p F takes 
the form: 

= C F [p F (t)} = -iu; F [^a,p F (t)]+L F [p F (t)] . (3.5) 

The study of the solutions of this equation in the holomorphic representation requires 
the introduction of Grassmann variables 6>, £, . . ., that anticommute with the operators 
and ct, and such that: 

0£ = -£0 , 9 2 = e = . (3.6) 

The elements \tfj) of the Hilbert space 7i F are now holomorphic functions ift(9) of the 
variable 9. However, since 9 2 — 0, their Taylor expansion contains only two terms: ifi(9) = 
^0 + ipi 9, with ifjQ and ipi complex parameters; they clearly represent the components of 
|^) along the vacuum and one-fermion states. ^ 

The inner product of two states \<p) and involves the integration over anticommut- 
ing variables (Berezin integral), defined by the conditions / 9d9 = 1 and J d9 = 0:[26] 

{^) = J V>* (0)0(0) e~ §e d9d9 , (3.7) 

where if}*(9) = ipQ + 9ip\ is by definition the adjoint of ip(9). 

Similarly, to an operator O acting on TC F there corresponds a kernel 0(9, 9); the result 
of its action on the vector \ip) is given by: 

m = J 0(9, ^(o e-ft de de . (3.8) 

Note that the identity operator is represented by the kernel e 9d . Furthermore, in this frame- 
work the fermionic creation and annihilation operators are realized by left multiplication 
and differentiation with respect to 9: 

d 

« f - 9 , a - , (3.9) 
so that «t is indeed the adjoint of a with respect to the inner product in (3.7). 



t Since the fermionic oscillator is a two-level system, a simple correspondence between 
the holomorphic and the standard matrix representation can easily be established; however, 
working with the holomorphic representation is in general more convenient, since explicit, 
closed expressions for p F can always be given, even in presence of n degrees of freedom. 
See also the discussion in Sect. 5 
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As in the bosonic case, since C F in (3.5) is quadratic in the operators (3.9), the kernel 
p F {9, 0\t) representing the state p F of the system can be taken to be of Gaussian form: 

p F {9,e-t)= 1 {t)e-~ er ^ . (3.10) 

For simplicity, also in this case we assume (a^) = (a) = for all times, so that terms 
linear in 9 and 9 are absent in (3.10). A more general Ansatz for p F (9,9;t) is discussed in 
the Appendix. The normalization condition: 

Ti[p F (t)] = J p F {9,9;t) e° e d9d9=l , (3.11) 

readily implies: 7 (£) = [l — -T(t)] \ so that the kernel p F in (3.10) contains only one 
independent unknown function. It can be conveniently recast in the following form: 

p F (9,9;t) = 1 (t)+[l- 1 (t)}99 , (3.12) 

explicitly showing that 7 and 1 — 7 represent the two eigenvalues of p F J Finally, the 
physical meaning of 7 (i) can easily be derived: 

0= 9 p F {9, 9; £)] e §9 d9 d9 = 7 (t) . (3.13) 

Insertion of (3.12) in the evolution equation (3.5) allows deriving the equation satisfied 
by the unknown function 7(t): 

7 (£) = -2 ( VF + a F ) 7 (t) + 2 VF , (3.14) 
whose general solution is simply: 

7 (t) = e -2(i*+<")*( 70 - 7oo ) + 7oo , (3.15) 
where 70 = 7(0) is the initial condition, while 

7oo = ■ (3.16) 

Since rj F and a F are positive constants, both 7 (i) and 1 — 7 (t) are non negative, so that 
< 7(t) < 1. Furthermore, independently from the initial condition, the density matrix 
p F describing the state of the fermionic oscillator always approaches for large times the 
equilibrium configuration: p|? = 7^ + (1 — 7oo)##; this is a thermal state, provided the 
inverse temperature f3 is introduced via the relation: e^ F = i] F /a F , with r\ F > a F . 



t Note that this simple rewriting of the Gaussian Ansatz is possible only in one dimen- 
sion; in presence of n degrees of freedom, the covariance r would be an n x n hermitian 
matrix and the Taylor expansion of (3.10) would be much more involved. 
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The evolution towards equilibrium is not in general associated with a monotonic in- 
crease of the von Neumann entropy S[pp] = Sp. Its explicit expression can be computed 
using the definition (2.26): 

S F (t) = - 7 (f) 1117(f) " [1 " l(t)} In [1 - 7(f)] , (3.17) 

while its time derivative reads: Sp = 7[hi(l — 7) — In 7]; one can easily check using (3.14) 
and (3.15) that Sp is always negative when 7 lays between 1/2 and 7^, while it is positive 
outside this interval. 

More precisely, as a function of 7, Sf grows from its minimum value Sf = at 7 = 
up to its maximum Sp = In 2 reached for 7 = 1/2, and then decreases, becoming again 
zero at 7 = 1. Therefore, Sp monotonically grows only when 7(f) increases in the interval 
[0, 1/2], or decreases in the interval [1/2, 1]. For rjp < <jp, this happens when 70 < 7^; 
indeed, this implies 7(f) < 700 < 1/2 and 7(f) > for all times. Similar conditions hold 
when rjp > o>; in this case to obtain a monotonic increase of entropy, one has to choose 
7o > 7oo so that 7(f) > 7^ > 1/2 and 7(f) < for all f. 

The case t\f = of is somehow special, since now S F > independently from the 
choice of the initial state; the density matrix p_p asymptotically approaches the infinite- 
temperature, totally disordered state p|? = e ee /2, for which the entropy is maximal, 
S F = In 2. 

As a final remark, note that in the case of the fermionic oscillator the sufficient con- 
dition for entropy increase discussed at the end of the previous section does not lead in 
general to useful constraints. Indeed, the inequality in (2.32) gives now the condition: 

S F {t) >2{ VF - <tf) [1-2 7(f)] . (3.18) 

Unless r\p = ap, the r.h.s. of this inequality becomes always negative for large enough 
times, as it can be easily realized by substituting for 7(f) its asymptotic value 7^. 



4. THE SUPERSYMMETRIC OSCILLATOR 

We shall now discuss the behaviour of an oscillator composed by both bosonic and 
fermionic degrees of freedom in interaction with an environment, under the hypothesis 
that its evolution is described by a quantum dynamical semigroup. The density matrix p 
representing the state of the system is now an operator on the Hilbert space H = Hb®'Hf- 
Its time evolution is described by an equation of the form (2.1), where both the total 
hamiltonian H and the dissipative piece L[p] are expressed in terms of bosonic, , a, and 
fermionic, cr , a, creation and annihilation operators, obeying 

[a f , cJ] = [a*, a] = [a, a 1 } = [a, a] = , (4.1) 

together with the standard commutation, anticommutation relations. 
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The hamiltonian H — Hb + Hp, the sum of the bosonic and fermionic terms of the 
form (2.9) and (3.2), possesses an additional property when the two frequencies are equal: 
ujb = ojp = oj. Indeed, the following charges: 

Q+ = u 1/2 aat, Q- = u 1/2 a ] a , (4.2) 

commute with the hamiltonian H = uj(a)a + a^a), and further: 

{Q+,Q_} = #, Q 2 + = Q 2 _=0. (4.3) 

This is the simplest example of a supersymmetry algebra. The system described by H is 
therefore supersymmetric and the conserved supercharges Q+ and Q- exchange bosons 
and fermions; further, from the algebra (4.3) one deduces that the ground state of H is a 
zero energy singlet and that all excited states form degenerate doublets. 

The additional piece L[p] in the evolution equation (2.1) will be taken to be the sum of 
the bosonic Ls[p] and fermionic Lp[p] linear operators already introduced in the previous 
sections. This is a natural choice since it assures integrability of the time evolution (L[p\ is 
again at most quadratic in the creation and annihilation operators) , while avoiding mixings 
between bosonic and fermionic degrees of freedom induced by the dissipative term; in other 
terms, L[p] is thus bosonic in character. 

Nevertheless, this simple form of L[p] does not in general assure supersymmetry invari- 
ance. In ordinary Quantum Mechanics, to an invariance of the hamiltonian there always 
correspond a conservation law and viceversa. For time evolution generated by equations of 
the form (2.1) this is usually not true: charge conservation and invariance (or symmetry) 
give rise to two different and in general unrelated conditions. 

To further elaborate on this point, notice that to the evolution equation (2.1) for the 
density matrix p there corresponds an analogous evolution for any operator X representing 
an observable of the system: 

^X = C*[X}=i[H,X} + L*[X] , (4.4) 

where the linear operator L* is the "dual" of L and it is defined via the following identity 

Tr(L*[X)p)=Tv(X L[p}) . (4.5) 

Consider now a symmetry of the hamiltonian H generated by the charge G, inducing the 
following transformation on the observables: 

X -> X' = U~ X XU , U = e iG . (4.6) 

This transformation will be an invariance of the system only when it is compatible with 
the evolution equation (4.4), i.e. C*[X'\ = U' 1 C*[X]U, for any X; equivalently, in 
infinitesimal form: 

[G,L*[X]]=L*[[G,X]} . (4.7) 
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This condition is clearly distinct from the relation that guarantees the time conservation 
of the mean value (G) = Tr[Gp(t)] of the generator G; recalling (2.1) and (4.5), from the 
condition d/dt(G) = for any state, one readily derives: 

L*[G] = . (4.8) 

In the case of the supersymmetric oscillator, the dual map L*[X] = + L^f-X"] 

can be easily deduced from (2.11) and (3.4). Explicitly, one finds: 

L* B [X] = tie (2 a^Xa - Xa^a - ataX) + a B (2 aXa j - Xaa j - aa) X) 

+ X B (2 a f AV - Xa ]2 - a ]2 X) + X* B (2 aXa - Xa 2 - a 2 X) , 

Lp[X] = r] F {2 a^Xa - I« f a - cfiaX) + o F (2 aXa 1 - Xaa j - aa^X) 
+ 2A f a t la t + 2X* F aXa . 



(4.9) 



The parameter of a supersymmetry transformation is anticommuting, so that the corre- 
sponding generator takes the form G = £ Q+, where £ is a Grassmann variable, commuting 
with bosonic operators, but anticommuting with the fermionic ones. Inserting it in (4.7) 
and using (4.9), after some algebraic manipulations one gets the following condition: 

(vb -(TB-VF + <tf) [X, G] + X* F {X, G^} = . (4.10) 

Since this relation must be true for any observable X, supersymmetry invariance is com- 
patible with the time evolution only when: 

Vb - &b = Vf - of , X F = . (4.11) 

The holomorphic representation is again particularly useful in order to discuss the 
behaviour of the state p(t) of the supersymmetric oscillator. The elements of the Hilbert 
space Ti will be now represented by holomorphic functions of the complex variable z and 
of the Grassmann symbol 9, while creation and annihilation operators will act on them 
following the rules in (2.16) and (3.9). The density matrix p will be now a kernel p(z, z; 9, 9), 
whose explicit expression can be taken to be of Gaussian form.t It can be expanded as: 

p(z, z- 6, 9) = p (z, z) + Pl (z, z) 99 . (4.12) 

The normalization condition Tr[p] = 1 now involves both ordinary and Grassmann inte- 
grals: 

J p(z,z;9,9) e~~ zz e dS dzdz d9 d9 = J [p (z, z) + Pl (z, z)\ e~~ zz dzdz = l . (4.13) 

Inserting the Ansatz (4.12) into the evolution equation for p allows deriving the fol- 
lowing conditions on the bosonic kernels P o and p±: 

p (t) = C B [po(t)} + 2 [t] F Pl (t) - o F po(t)] , (4.14a) 
p 1 (t) = jC B [p 1 (t)]-2[rj F p 1 (t)-a F p (t)] , (4.146) 



t Here again we assume vanishing initial averages (af), (a), (a^), (a) 
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where the linear operator £b[p] is as in (2.10). It follows that the combination p + p\ 
satisfies the same evolution equation discussed in Sect. 2 for the case of a single bosonic 
oscillator. The Gaussian Ansatz p B (z,z) in (2.17) can equally well be adopted here for 
Po + Pi) since performing the Grassmann integrations, one consistently finds (compare with 
(2.20)): 

/q2 
— [p (z,z;t)+p 1 (z,z;t)] e~~ zz dzdz = x{t) , 

(a t2 )(t) = Tr[a t2 p(t)] = J z 2 [p (z, z;t) + p^z, z;t)] e~~ zz dz dz = x(t) , (4.15) 

—z[po(z,z;t)+ Pl (z,z;t)] e~ zz dzdz = y(t) . 

As a consequence, the time evolution of these quantities is that given in (2.22) and (2.24). 
Inserting back this result into (4.14a), one obtains: 

Po(t) = C B [p (t)} -2( VF + <r F )po(t) + 2 VF p B (t) . (4.16) 

The form of this equation suggests to look for a solution in which po(t) differs from psif) 
by an unknown multiplicative function 7f(^)- It can be identified with the function 7(f) 
studied in the previous section, since it satisfies the same equation (3.14) and has the same 
physical meaning: 

(<W)(t) = Tr [aa* p(t)] = J p (z,z;t) e~ zz dzdz = ^ F {t) . (4.17) 
As a consequence, p\ = (1 — 7f)p_b 5 and therefore one finally finds: 

p(z,z;6,6) = [ 7 f + (1-7f)^] Pb(z, z) = p F (6, 9) p B (z,z) . (4.18) 

Not surprisingly, the density matrix that solves the evolution equation (2.1) in the case 
of the supersymmetric oscillator is in factorized form; its behaviour can be deduced from 
the analysis of the previous sections, provided the conditions (4.11) for supersymmetry 
invariance are taken into account. 

In particular, p approaches an equilibrium state for large times only when t\b > ct b , 
which also implies: r\ F > o F . This limiting state is thermal, with inverse temperature (3, 
only for Xb = and t^b/^b = Vf/cf = , which implies, recalling the condition (4.11): 
Vb = Vf and ob = cr F .[30] 

Also in the case of the supersymmetric oscillator, the total entropy does not have 
in general a monotonic behaviour during the approach to equilibrium. Since the density 
matrix p is in factorized form, the total entropy S will be the sum of the bosonic and 
fermionic contributions. Using the variable 75 = % — 1/2 > 0, where x is defined as in 
(2.30), and recalling the results of the previous sections, one has: 

S = (75 + l) In (75 + l) - 7s In 7s - (7^ - l) In (7^ - l) - 7^ In 7^ . (4.19) 
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Its time derivative, that can be expressed as: 



S = 7B ln^l + ^+ 7F ln^l--^ , (4.20) 

does not have in general a definite sign, although possible compensations between the 
bosonic and fermionic contributions can concur to a positive r.h.s. for certain time inter- 
vals. 

As discussed at the end of Sect. 2, a bound on S can be obtained by working directly 
with the definition (2.26) and the equation (2.1). In the present case, this procedure gives: 

S> 2(a B - VB ) +2(a F -rj F )[2 lF -l] = 4(a F - Vf ) lF , (4.21) 

where the identity is a consequence of the condition (4.11). Since < 7f(£) < 1, the 
inequality (4.21) assures S > for a F > r\ F . However, this condition would lead to a 
rather singular behaviour for the bosonic part of the density matrix in (4.18), and thus 
for the whole p. In fact, also o~b would be greater than rjs and, as discussed in Sect. 2, 
this implies an infinitely growing average occupation number. In conclusion, although 
inducing a partial compensation between the bosonic and fermionic contributions to S, the 
supersymmetry condition (4.11) is in general not enough to guarantee monotonic entropy 
increase for all times during the evolution of the system. 



5. DISCUSSION 

All the considerations developed in the previous sections for single oscillators can be 
generalized to the case of n independent oscillators, both bosonic and fermionic. Their 
interaction with an external environment can still be consistently described in terms of 
quantum dynamical semigroups, so that their time evolution can be modelled by means 
of equations of the form (2.1), (2.3), with operators Lk linear in the relevant fundamental 
variables. However, the coefficients r and s of (2.4) become now matrices, and the number 
of independent constants characterizing the dissipative part L[p] rapidly increase with n, 
making the evolution equation (2.1) rather involved. 

Nevertheless, various simplifying conditions can be imposed to reduce, at least in part 
this arbitrariness. Those involving symmetry properties are the most physically interesting. 
As discussed in the Introduction, the interaction between system and environment can 
be considered in general to be weak; therefore, in many instances, the presence of the 
environment should not to be able to alter the symmetry properties of the system. In the 
case of n isotropic oscillators, the hamiltonian is invariant under the action of the group 
SU(n); it is then quite natural to assume the same invariance property to be valid for 
the full evolution equation. As discussed in Sect. 4, this can be achieved by imposing the 
condition (4.7) for any element G of the SU(n) algebra. 

Also in this more general setting, the holomorphic representation appears to be a 
particularly convenient framework to analyze the behaviour of the solutions of (2.1). It 
requires the introduction of n commuting or anticommuting complex symbols, that allow 
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realizing the corresponding creation and annihilation operators as multiplication and dif- 
ferentiation by these variables. The kernel representing the system density matrix can still 
be taken to have a generic Gaussian functional expression. However, the various "coeffi- 
cients" in the exponent, suitable generalizations of the functions x, x and y in the bosonic 
case and of r in the fermionic one, become now n x n matrices. They obey quadratic 
(Riccati-like) time evolution equations, whose solutions can always be obtained, albeit in 
general in terms of implicitly defined functions. [31, 32] 

Although developed in the analysis of simple open systems, the techniques described 
in the previous sections are actually very general; they can be used to study the dynamics 
of more complicated models, for which the operator C in (2.1) is not quadratic in the 
relevant variables. In these cases, complete explicit expressions for the density matrix p 
as solution of (2.1) can not in general be given. Nevertheless, approximate expressions for 
p, typically in Gaussian form, can be obtained via the application of suitable variational 
procedures. 

Indeed, equations of the form (2.1) can be derived by mean of a suitable variational 
principle, [33] obtained by generalizing the one yielding the Liouville - von Neumann equa- 
tion in ordinary Quantum Mechanics. [28, 34] In the case of isoentropic time evolutions, 
these variational techniques have allowed detailed discussions of a wide range of physical 
phenomena, from statistical physics to inflationary cosmology. [34, 35] Their application to 
the study of quantum dynamical semigroups within the framework presented in the previ- 
ous sections will surely provide new insights on the behaviour of open quantum systems. 



APPENDIX 

The Gaussian kernels pB(z,z;t) and pp(0,9;t) representing the density matrices for 
the bosonic and fermionic oscillators discussed in Sect. 2 and 3 lead to vanishing averages 
for the corresponding creation and annihilation operators. This condition can easily be 
released by introducing a more general Ansatz. 

In the bosonic case, take: 

p B (z, Z] t) = 1 c -W77t[ 2 ^) [z-v{t)][z-v(t)]-x(t) [z-v(t)] 2 -x(t)[z-v(t)] 2 ] +zz ^ 

that differ from the expression (2.17) because of the presence of the two additional functions 
v(t) and v(t). Trace conservation for all times, Tr[p B (t)] = 1, still implies: 

N(t)=y 2 (t)-\x(t)\ 2 , (A.2) 

while hermiticity requires: v(t) = With this choice for p B (z,z;t), the averages of 

and a are in general non vanishing: 



(a f )(t) = Tr[a f p B (tj] = J zp B {z,z;t) e~ Sz dzdz = v(t) , 

— [p B (z,z;t)] e~~ zz dzdz = v(t) . 



(A3) 
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The time evolution equation (2.10) implies the following homogeneous equation for v(t) 

v(t) = 2(7 lB -a B + ioJ B ) v(t) , (AA) 
so that v(t) is non vanishing only if its initial value v(0) is different from zero: 

v (t) = e -(vB-<r B +iu B )t v ^ _ 

The physical meaning of the remaining functions x, x and y appearing in (A A) is 
slightly changed with respect to those studied in Sect. 2; now one finds: 

x = (a 2 ) - (a) 2 , 

x = (a t2 > - (a f > 2 , (A6) 
V = (aa 1 ) - (a) (a f ) . 

Nevertheless, one can check that these functions still obey the evolution equations (2.21), 
so that the considerations and the discussions of Sect. 2 apply to this more general situation 
as well. 

In the case of the fermionic oscillator, the most general Gaussian Ansatz for the kernel 
Pf(0, 0', t) can be written as: 

MM;*) = 7(f) e-J* [**<*>»-*<*) . (A7) 

The normalization condition Tr [/9p(t)] = 1 gives A(t) = ^(t) — 1, while hermiticity implies 

<p{t) = Mt)}*. 

By performing the integration over the anticommuting variables, one finds that the 
function ^(t) retains its meaning as (a a*) also in this more general setting, and therefore 
still obeys the evolution equation (3.14). 

On the other hand, the two additional functions (p(t) and <p(t) in (A. 7) represent the 
averages of and a, 



(a f )(t) =Tr[cJp F (t)] = J 8p F (8,e;t) e Se d6d9 = p(t) , 

(a)(t) = Tr[ap F (t)] = J £[p F (W)] e §e d9 d9 = <p(t) 

and, as a consequence of (3.5), obey the following evolution equations: 

(p(t) = -2 (t]f + cr F + ioj F ) <p(t) + 2X F <p(t) , 
(p{t) = -2 (r] F + a F - iu F ) (f(t) + 2X* F (f{t) . 

The general solution is given by: 



(A.8) 



(A.9) 



cos(Q F t) — — — sin(Q F t) 

{It? 



V (0) + ^sm(fl F t)m} , (AAO) 



(p(t) = e - {r > F+aF)t 

1 /2 

where fl F = [u F —4\X F \ 2 ] for u F > \X F \. Hyperbolic functions appear in the expression 
(A. 10) when u F < \X F \; however, thanks to the inequality |A^| 2 < r\ F a F (compare with 
(2.13)), (f(t) always vanishes for large times. 
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